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It is estimated that half the world's population cooks over an open biomass fire; improved biomass cooking 
stove programs have the potential to impact indoor air quality, deforestation, climate change, and quality of 
life on a global scale. The majority of these cooking stoves operate in a natural convection mode (being driven 
by chimney effect buoyant fluid forces). Simplified theories for understanding the behavior of this 
unexpectedly complex combustion system, along with practical engineering tools to inform its design, are 
markedly lacking. A simplified model of the fundamental stove physics is developed for predicting bulk flow 
rate, temperature, and excess air ratio based on stove geometry (chimney height, chimney area, viscous and 
heat release losses) and the firepower (as established by the stove operator). These parameters are intended 
to be fundamental inputs for future work understanding and improving biomass cook stove emissions and 
heat transfer. Experimental validation is performed and the simplified model is shown to be accurate and 
applicable to typical stove operation. Carbon monoxide and particulate matter emissions data have been 
recorded in conjunction with the validation data. The initial results are presented and indicate that the excess 
air ratio is a promising tool for reducing carbon monoxide emissions. 

© 2011 International Energy Initiative. Published by Elsevier Inc. All rights reserved. 


Introduction 

It is estimated that about half the world's population—more than 
3 billion people—cook over an open biomass fire. This energy source is 
especially prevalent in rural areas of developing countries, where an 
estimated 90% of households rely on biomass fuels for cooking (de 
Koning et ah, 1985). Despite being a simple, convenient, and 
potentially carbon neutral energy source, several important chal¬ 
lenges are apparent. 

Combustion of solid biomass fuel is especially prone to incomplete 
combustion, often resulting in harmful emissions. These emissions, 
when combined with poor ventilation, often lead to poor indoor air 
quality and a number of adverse health effects. For a more complete 
description see de Koning et ah (1985); Naeher et al. (2007) and 
Smith et al. (2003). Over 1.6 million deaths a year are attributed to 
poor indoor air quality as a result of domestic biomass combustion 
products (Smith et ah, 2003). Poor indoor air quality is thought to be 
the leading cause of death worldwide for children under the age of 
five (Bruce et ah, 2000). 

Beyond the direct exposure health concerns—and with such 
extensive use—cooking with biomass fuel can also have important 
global implications on greenhouse gas and black carbon emissions. 
Smith (1994) and Haines et al. (2010) provide a good discussion of the 
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overall global warming potential; cooking with biomass fuels is 
thought to account for 1-3% of all human generated global warming. 
In Venkataraman et al. (2005), biomass cooking is suggested to be the 
dominant source of global, energy related, black carbon emissions. 
Biomass cooking is also a significant source of deforestation and 
consequent climate forcing as discussed in Mwampamba (2007). For 
many of the world's poorest groups, acquiring biomass fuel also 
presents a heavy burden of time and/or money. 

Biomass (improved) cooking stoves 

Biomass based cook stove technology is increasingly being used to 
address these issues (Barnes et al., 1993; Baldwin et al., 1985). Cook 
stove performance is typically characterized by the following two 
metrics: reduction of harmful emissions, primarily carbon monoxide 
(CO) and particulate matter (PM), and increased thermal efficiency. 
These are not the only concerns, but are readily quantifiable metrics 
used often by the stove designer. 

The majority of these biomass cooking stoves are wood burning 
and operate in a natural convection mode—being driven by chimney 
effect buoyant fluid forces. This configuration will be the focus of this 
research, however the results are likely applicable to other stove/fuel 
configurations. 

The “three stone fire” is the most basic, but still extremely 
common, method for cooking with biomass fuel, and is often used as a 
benchmark for comparison. The three stone fire is simply an open fire 
—gaining its name from the stones used to hold thecook piece over the 
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fire. Biomass cook stoves have been able to reduce CO and PM 
emissions by 50-75%, and increase fuel efficiency by 30-50%, 
compared to the three stone fire configuration (Jetter & Kariher, 
2009). These reductions are accomplished largely by enclosing the fire 
in a combustion chamber. The primary benefits of chamber enclosure 
are airflow control (providing better mixing with reduced quenching) 
and reduced radiative and convective heat transfer losses. The 
combustion chamber also allows the stove designer to locate the 
cook piece for increased thermal efficiency. 

There exists a wide variety of stove designs with varying 
performance (Smith, 1989; MacCarty et al., 2010). Stove efficiency 
and emissions are very sensitive to the combustion chamber shape, 
material, chimney height, chimney diameter, and cook piece 
placement. These design choices have been studied by “improved” 
cook stove projects for over 40 years. A great deal has been learned, 
resulting in empirically/observationally established design guidelines, 
but rarely in theory based analytical or numerical design tools (Smith, 
1989). The well known “rocket elbow" stove design, shown in Fig. 1 is 
a good example. The design has been empirically studied and refined 
for at least two decades, and the resulting design guidelines have been 
documented (Baldwin, 1987; Hudelson et al., No Date; Scott, 2003; 
Bryden et al., 2005). Moving beyond this kind of empirical/ 
observational approach, the complexity of the physical and chemical 
processes, as well as their interconnection, increase rapidly. 

For the development of analytical theory and practical engineering 
tools, it becomes necessary to extract the simplified fundamental 
processes governing basic stove behavior. In the field of biomass 
cooking stoves research, this type of simplified, foundational theory is 
notably lacking. 

A simplified stove flow model 

A simplified model of the fundamental stove physics will be 
developed for predicting bulk flow rate, temperature, and excess air 
ratio, based on stove geometry (including chimney height, chimney 
area, viscous and heat release losses) and the operating firepower. 
Stove geometry inputs are of interest to the stove designer as they are 
basic tools for influencing stove performance. The stove firepower is 
primarily established by the operator, and is consequently also 
included as a model input. A given stove/fuel configuration can only 
be operated within a certain firepower window. Within this window 



the stove operator establishes the operating firepower of the stove by 
controlling the fuel feed rate and positioning. Excess air ratio will be 
shown to provide some insight into the upper limit of this range. 

The resulting bulk flow rate, temperature, and excess air ratio 
behavior are intended to provide insight and fundamental inputs for 
beginning to understand, predict, and improve stove performance. 
Bulk flow rate and temperature are important inputs for analytical, 
empirical, or numerical studies of the thermal efficiency and 
impinging flow heat transfer to the cook piece (pot). A simplified 
model is also conducive to use in optimization and automated 
routines. Excess air ratio and the oxidizer co-flow rate affect 
combustion chemistry and aerodynamics, and thus the subsequent 
emissions and heat transfer. Although excess air ratio is most 
traditionally applied to premixed combustion, it has also been applied 
to biomass boiler design and operation for affecting emissions. Initial 
results included here will suggest its suitability for biomass cookstove 
design. 

Model setup 

A schematic setup of the analytical stove flow model to be 
developed is shown in Fig. 2. The 1-D flow shown in Fig. 2 is assumed. 
The physics of the stove have been simplified to two basic and 
fundamental driving processes: heat addition from combustion (at 
point 2), and kinetic energy addition/conversion (between points 1 
and 2) due to the chimney effect. The two processes are inter¬ 
connected and together govern the simplified overall stove operation. 

Heat addition 

The driving fire releases heat into the passing stove flow, 
increasing its temperature and decreasing its density. Heat addition 
is assumed to be perfectly efficient and is greatly simplified as taking 
place entirely and instantaneously at point 2. For a given mass flow 
rate, m A , and heat release rate (essentially the current stove 
firepower), Q in , the bulk flow temperature increase can be calculated 
using the first law of thermodynamics for a control volume. For an 
isobaric system with no mechanical work, assuming ideal gas 
behavior, and constant potential and kinetic energy, the temperature 
increase can be calculated using the constant pressure specific heat of 
air (c p ). The heat addition from combustion translates into an 
enthalpy increase ( h c to h H ) distributed over the mass flow rate 
( m A ) as shown in Eq. (1). 

Qin = m A (h c -h H ) 

= m A j TH c(T)dT (1) 

J Tc 

^A^p,avg(^H T c ) 



Fig. 1. Insulated rocket elbow stove. 


Fig. 2. Simplified model setup. 
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Where T is the bulk flow temperature, and the subscripts H and C 
denote hot and cold respectively. A constant specific heat, evaluated at an 
average temperature, will be assumed. Heating takes place over a range 
from approximately 300-1300 K, over which specific heat varies by a 
maximum of 15%. A temperature dependent specific heat can be 
calculated using a piecewise polynomial of the form described in (Kee 
et al., 1990), the resulting error never exceeds 2%, and is considered 
negligible. 

Chimney effect 

Airflow through the stove is dependent on the chimney effect- 
resulting from the buoyant force of the decreased density air after 
combustion heat addition. The chimney effect equation will be briefly 
derived. The decreased density air in the chimney between points (2) 
and (3) creates a lowered pressure at point (2) compared to the 
ambient pressure at point (1). In a fluid, pressure increases with 
depth, and the pressure increase (AP) can be calculated using the 
hydrostatic equation, given below in Eq. (2). 


difference using the integral form of Bernoulli’s equation for 
compressible flow, as shown in Eq. (5). 

A Pi -2 = gh(pAmb~PH ) 

_ 1 2 ( 5 ) 

— 2ptf V 2 


Where the hot gas density and velocity (v) at point (2) are used for 
calculating the stove flow kinetic energy. Velocity, volume (V) and mass 
flow are related by the flow area (A) and density. Again assuming ideal gas 
behavior, the density is related to the temperature by the ideal gas law. 
Utilizing these relations and solving Eq. (5) for volume and mass flow rate 
yields Eqs. (6)-(9). Eq. (8) is the commonly known form of the chimney 
effect equation. The mass flow rate form shown in Eq. (9), will be a more 
applicable and interesting form for the purposes of this paper. 


(PAmb-pH)gh 




A P = gjp(h)dh 





P Amb Ph 
P H 



Where g is gravity, and the density (p) is a function of chimney 
height (ft). The pressure at points (1) and (2) can be calculated using 
Eq. (2), the former by starting at point (3), and integrating while 
traveling down the chimney to point (2), and the latter by traveling 
outside of the chimney through ambient air (with constant ambient 
density (p Am b) to point (1). The difference is found by subtracting the 
two results, as shown in Eq. (3). 

AP,_ 2 = (p 3 + g J1 PAmbdh') 

~(P 3 +gj!p(h)dh) (3) 

= gPAmb h ~gjlP( h ) dh 

If we assume the chimney walls to be adiabatic, than the 
temperature and density of the gas in the chimney (T H and p H ) will 
remain constant, and Eq. (3) simplifies to Eq. (4). 

APi-2 — gh{PAmb~ Ph) (4) 

In reality, approximately one third of the energy may be lost to the 
chimney walls. Including this effect, the integral from points (2) to (3) 
could be performed if the temperature/density profiles were known. For 
the simple model presented here, no attempt is made to calculate these 
profiles, and the adiabatic chimney assumption is maintained. 

Chimney flow results from the pressure difference calculated in 
Eq. (4); the gain in kinetic energy of the chimney flow at point (2) 
from the stagnant ambient air at point (1) is equated to this pressure 







Where the loss coefficient (C), has been introduced to account for 
uncertainties and inefficiencies in the chimney effect—primarily including 
viscous losses, chimney wall heat transfer, and the unrealistic ideal point 
heat addition at point two. With regards to the point heat addition 
simplification: In reality, combustion heat is released throughout the 
chimney at the transiently located flame front; when this heat addition 
takes place after point 2, the chimney density profile used in Eqs. (3) and 
(4) is less effective in producing the chimney effect in Eqs. (8) and (9). In 
this case, as well as with viscous and chimney wall heat losses, the result is 
a reduced chimney effect and 0 < C<1. Any or all of these losses may be a 
function of flow rate or temperature but for now this dependence is 
unknown and neglected. 

The resulting two equation system 

Solving Eqs. (1) and (9) simultaneously, a simple model for 
predicting stove bulk flow and temperature is determined. The system 
of equations is shown in Fig. 3, along with the input, output, and 



Results in... 
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• Mass flow rate, m.\ 

• Air Fuel Ratio or Excess Air Ratio 
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• Thermal Efficiency 

• Flame Shape 


Fig. 3. Overall system of equations with inputs and outputs. 
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feedback connection between the combustion heat addition and the 
chimney effect flow. Stove geometry and firepower are inputs as 
previously discussed. This system will be further developed into a 
dimensionless form in Part 2, Agenbroad et al. (2011), with 
advantages including scale similarity and a reduction of experimen¬ 
tally independent parameters. 

The two equations can be solved iteratively, in this case using a 
MATLAB® script. The iterative solution results in a steady state bulk 
flow rate and temperature for a given stove geometry and operating 
firepower/fuel feed rate. Although the actual stove behavior is 
inherently transient, these steady state values will be shown to 
have physical significance for temporally averaged behavior. 

The air fuel ratio (AFR) and excess air ratio (EAR) are calculated 
using the chimney effect driven mass flow rate and the fuel mass burn 
rate. The mass flow rate of fuel (rh F ) is calculated from the firepower 
(Q.,n) and the heating value (HV) of the fuel as shown in Eq. (10). The 


AFR and EAR can then be calculated as shown in Eqs. (11 )— (13), 
where AFR stoich is the AFR for stoichiometric combustion. 

• _ Qin 

mp ~ LHV 

(10) 

AFR = 

m F 

(11) 

r jh _ AFRstoich 

AFR 

(12) 

%EAR = (->-« 

0 

(13) 

The lower heating value (LHV) 

is used—where the LHV is 


essentially the heat of combustion. The higher heating value is not 
be used because the latent heat of the water vapor is significant, and 
will not be recovered from the flow. 


Combustion, heat transfer, and viscous losses were neglected and the 
loss coefficient was not used (LC = 1). Volume flow rate and 
temperature are as expected; increasing firepower leads to increasing 
flow rate and temperature. Mass flow rate behavior is more interesting; 
a peak mass flow rate is observed at a relatively low firepower, after 
which mass flow rate decreases with increasing firepower. The 
increasing volumetric flow rate is essentially competing with the 
decreasing flow density due to the increasing flow temperature. 

Air/Fuel and excess air ratio 

As firepower is increased, this peak, and the following reduction, in 
mass flow rate have an interesting effect on excess air ratio behavior. 
Excess air ratio (as calculated from Eq. (13)) is shown with temperature 
and mass flow rate behavior in Fig. 5. The excess air ratio is large while 
operating at low firepowers, but decreases rapidly as firepower is 
increased. 

The x-axis range shown in Fig. 5 spans the practical operating range 
for most typical single family sized cooking stoves. At high firepower, 
relatively low excess air ratios are approached (although 100% excess air 
ratio would still be considered high for boiler applications). A theoretical 
upper limit for firepower (the stove's maximum firepower) can be 
predicted at 0% excess air ratio. At this point exactly sufficient oxidizer is 
supplied for combustion of the current fuel supply rate. A further 
increase (assuming no limit in pyrolysis fuel supply is reached first) will 
not have sufficient oxygen for combustion, and a corresponding increase 
in stove firepower output cannot be achieved. However, considering the 
poor mixing of actual stove operation, this limit will be reached at some 
point before 0% excess air ratio. 

It is unlikely that high excess air ratios should limit a minimum 
firepower. For the case of the “three stone fire,” excess air ratios would 
be considered to be very large, if not infinite, for the entire range of 
operation. High excess air ratios may however lead to increased 
quenching and cooler combustion temperatures, likely promoting 
incomplete combustion and a reduction in thermal efficiency. 


Model results and discussion 


Geometric effects 


The system of equations described above has been solved iteratively 
over a typical range of operating firepowers and the results are shown in 
Figs. 4-6. Rocket elbow stove geometry used is given in Table 1 . A LHV of 
18.280 MJ/kg and a C 4A H 63 0 2 .5 fuel composition (ultimate analysis 
taken from Tillman (1979) for Douglas Fir) were used. 

Bulk flow rate and temperature 

Bulk flow temperature, volumetric and mass flow rate behavior are 
shown for constant geometry and varying firepower in Figs. 4 and 5. 


And finally, we look at the effect of varying geometry inputs. The 
three geometry inputs appear together in Eq. (9); their effect can be 
summarized by their product CAVh. The effect of doubling and 
halving the product of this group is shown in Fig. 6. Increasing the 
geometry product (either by increasing stove size or reducing losses), 
will increase both the mass flow rate and Excess Air Ratio (EAR), while 
decrease the Temperature (T), for a given firepower. Decreasing 
product has the opposite effect. 

These geometric inputs represent important design choices/tools 
to the stove designer. The stove designer typically cannot control the 




Fig. 4. Volume flow rate, mass flow rate, and temperature behavior. Fig. 5. Mass flow rate, temperature and excess air ratio behavior. 
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Fig. 6. The effect of varying stove geometry inputs. 
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users choice of firepower within the stove's operating range. It is 
therefore interesting to examine the effect of these geometric 
parameters over a range of stove operation as shown. 


Experimental validation setup 

Experimental validation has been performed using two sizes of the 
standard rocket elbow geometry fed by hand with wooden stick fuel 
for temporally averaged data sampling. 

In contrast to more strictly controlled/idealized operation, the 
model is validated for conditions similar to actual use—assessing its 
practical applicability, but requiring temporal averaging due to its 
inherently transient nature. A “test sweep” is performed consisting of 
a sweep through five target firepowers. Firepower is held as constant 
as possible for each of the five, 15-minute long, sampling periods. 
These 15-minute sampling periods are then temporally averaged for 
each measured parameter. Firepower is controlled by varying fuel 
feed rate, number of sticks, and stick positioning. A 15-minute 
transient period is allowed between sampling periods for the stove to 
adjust to each new firepower. The rocket elbow's low thermal mass 
and the fiberfrax insulation also help minimize transient effects. An 
example of the raw data collected for one test sweeps is shown in 
Fig. 7, with the warm-up, sampling, and transient periods shown. 

Two stove sizes (of the dimensions given in Table 1 ) were used to 
test the model's accuracy for varying geometry. Two test sweeps have 
been performed with each stove geometry. 

Validation was performed for somewhat idealized stove operation 
without a cook piece (pot) in place. This decision was made in an 
attempt to isolate basic simplified model physics—assessing baseline 
model accuracy while minimizing the effect of the unknown loss 
coefficient and neglected heat transfer. Part 2, Agenbroad et al. 
(2011), will extend and examine model application to the more 
practically applicable stove operation with a cook piece in place. 

By measuring stack exhaust temperature and %0 2 along with mass 
based stove emissions of CO and C0 2 , the firepower, mass flow rate 
and temperature are determined for comparison to the model 
predicted behavior. 


Measuring firepower 

Firepower is determined using a simple carbon balance of the 
measured CO and C0 2 mass emissions. Stove emissions are captured in 
a hood after 22 times volumetric dilution. Mass based CO and C0 2 
emissions are then measured using a Non-Dispersive Infrared (NDIR) 
Sensor. For the one step combustion reaction shown in Eq. (14) and 
the previously assumed fuel chemistry and heating value (LHV), the 
firepower is calculated as shown in Eq. (15). 


aC x H Y O z + b(0 2 + 3.76N 2 ) 
-^dC0 2 + eH 2 0 +/0 2 
+b(3.76N 2 ) + gCO 




^LHV in 



X 


^Molecular Weight in 


l<g \ 
kgmol J 


(14) 


(15) 


The single step reaction shown in Eq. (14) includes the combustion 
of gaseous volatiles as well as the surface oxidation of leftover carbon 
char. The two processes require different amounts of oxygen and have 
different heat release values. Some error results as the proportion of 
carbon originating from each process is unknown. However, this 
effect is minimized as wood volatiles contributions (on a mass basis) 
are much larger than those from char oxidation. Moreover, after 
temporal averaging their proportion remains relatively constant. 
Remaining error is neglected relative to other effects. 

Temperature 

Heitor and Moreira (1993) provide a discussion of in-flame 
temperature and species sampling. Bulk flow temperature is mea¬ 
sured using a K-type thermocouple placed at the approximate center 
of the chimney about 1 cm above the chimney exit. The technique 
works well for the temperatures range, location, and accuracy desired 
for this paper. 

Mass flow rate and %0 2 


Table 1 

Rocket elbow (RE) stove geometry parameters. 



4 in RE 

5 in RE 

Total height 

250 mm 

320 mm 

Chimney height 

200 mm 

255 mm 

Chimney diameter 

100 mm 

130 mm 


Stack exhaust volumetric %0 2 measurements are used in deter¬ 
mining actual stove mass flow rate. Mass flow rate is calculated from 
stack exhaust %0 2 as shown in Eqs. (16)-(18).Z0 2 measurements are 
taken using the 0 2 analyzer function of a N O x 5210 made by ECM (Los 
Altos, CA). The accuracy of the 0 2 analyzer was verified with a 3 point 
calibration prior to use. The sampling location is the same as the 
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Fig. 7. Example raw data before temporal averaging from 4 in rocket elbow Test 1. 


thermocouple and care is taken to sample at a rate no higher than 1% 
of total stove airflow to avoid affecting stove behavior. The analyzer 
measures wet %0 2 values, the dry values were arrived at by first 
passing the sampled airflow through a desiccant. 

Stack exhaust %0 2 is related to the single step reaction given in 
Eq. (14) as shown in Eq. (16). For a dry^0 2 measurement, e is equal to 
zero. By substituting element conservation equations taken from the 
reaction in Eq. (14) into Eq. (16), the molar rate of oxygen (b) can be 
calculated as a function of only the molar rate of fuel (a, calculated 
using Eq. (15)) and the fuel composition(C x HyO z ) as shown in 
Eq. (17). Finally, the mass flow rate of air is arrived at by multiplying 
by the molecular weight as shown in Eq. (18). 


point to be 13/im). Downstream of the cyclone, Whatman PTFE (Fisher 
Scientific #05 — 717 — 5) Teflon filters of 47 mm diameter and a 2/im 
pore size were used to collect PM. PM samples were taken after 
dilution at the top of the fume hood near the NDIR sampling point. 
Sampling area was about 1/680 of the total flow area. From previous 
experience with this setup, the limit of detection (LOD) has been 
calculated to be 15/ig, with a limit of quantification for a sample to 
have statistical significance at 51/ig. For this setup, PM measurements 
were greater than the LOD/LOQ, however some samples, if taken for 
less than 15 minutes would likely not reach the LOD/LDQ. 

Validation results and discussion 



/• 100 

d + e + / -T 3.76 b 



a(z- 

%0 2 -4.76—100 
x^50 -^-%0 2 -ax^j 


(16) 


(17) 


Experimental data taken from the four validation test sweeps are 
compared with model predictions in Figs. 8-12. Model behavior is 
shown for a loss coefficient (C) of 0.5 and for the ideal case with a loss 
coefficient of one. The two rocket elbow (RE) stove geometries are 
summarized Table 1. 

Mass flow rate and temperature 



k kgmol 
s 


(32 + 3.76-28) 


kgmol 


(18) 


Synchronizing instrumentation 

Sampling with the thermocouple is essentially instantaneous. 
Time delay for NDIR sampling with respect to the thermocouple was 
determined, through spiking, to be approximately 70s. The 0 2 sensor 
delay is estimated to be about 45 s, determined by comparison of the 
raw data between the thermocouple, NDIR, and the 0 2 sensor. 

Particulate matter 

Particulate matter data were measured for each 15-minute sample 
using gravimetric analysis. PM 10 data were taken by first using a URG 
Corporation 2000-30ENB cyclone to remove particles as near above 
the 10/im cut point as possible (specifications indicate the actual cut 


Fig. 8 shows the averaged %0 2 values resulting from the firepower 
sample points taken in each of the four tests with comparison to the 
model predictions. The corresponding mass flow rate values as 
calculated from the %0 2 measurements using Eqs. (17) and (18) are 
shown in Fig. 8. Error bars indicate one standard deviation in the 
firepower and %0 2 or mass flow values calculated from the 15-minute 
data sample. The experimental validation data agrees well with the 
model predicted behavior. 

The validation data shown in Fig. 8 are in agreement with the 
interesting model prediction of a decreasing mass flow rate with 
increasing firepower. Higher uncertainty is observed at lower 
firepowers. These firepowers are harder to create and maintain 
steadily. These lower firepowers also have decreased mixing and 
uniformity throughout the chimney. 

Corresponding validation temperature data are shown in Fig. 8. 
Trends are predicted well, but in this case the sampled data lies closer 
to the ideal model predictions. The lower than expected measured 
temperature is likely due to chimney wall radiative and convective 
heat loss. The loss coefficient used in Eqs. (8) and (9) only accounts for 
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• 4in. RE Test 1-4in. Model w/ LC 

o 4in. RE Test 2 - 5in. Model w/LC 
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Fig. 8. Temperature samples with comparison to model. 
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Fig. 9. Scatter plot of all %0 2 and temperature data. 


transient nature of the stove (evident for example in the raw data 
shown in Fig. 7). These plots provide further assurance that the model 
physics are applicable and raise the question of whether shorter, or no 
temporal average could be used. 

Stove firepower range 

The stoves were operated in a typical manner over the range of 
firepowers which seemed appropriate for cooking use. High firepower 
data points were determined by the highest practical fueling rate that 
could be conveniently maintained. The operating range of the 5 in. 
rocket elbow geometry is observed to be higher than that of the 4 in. 
geometry in all of the validation sampling. The larger stove has higher 
excess air ratios for a given firepower, and it is interesting to note the 
upper firepower limit occurred at very similar excess air ratios. For 


lowered chimney effect; the temperature is still calculated for the 
ideal case. A second loss coefficient can be introduced into Eq. (1) 
accounting for lost heat in the heat addition equation. Up to one third 
of the firepower of the stove is thought to be lost to the stove walls; if 
this heat loss equation coefficient is introduced and set to 0.65-0.75 
the model predictions agree well with the validation measurements. 

Instantaneous/un-averaged results 

The scatter plots shown in Fig. 9 were created using the raw, un¬ 
averaged, temperature, %0 2 , and corresponding firepower data. The 
plots include data sampled at 0.7 Hz from all four test runs, including 
transient periods between 15 minute data samples. The data have 
been corrected for the time delay between sensors as described under 
validation setup. The data fits well considering the inherently 



Fig. 10. Excess air ratio samples with comparison to model. 
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Fig. 11. Combustion efficiency averages for validation samples. 

ideal fuel/air mixing, the maximum stove firepower is predicted to 
occur at 0% excess air ratio. With the validation results, the maximum 
firepower was consistently observed at about 100% excess air ratio. 
After accounting for less than ideal mixing, this excess air ratio 
appears to inform an upper limit of a stoves firepower range. 
Pyrolysis, being endothermic, relies largely on radiative feedback 
from the flaming region. Prior to experimentation, it was considered 
whether pyrolysis rate might also limit the firepower range, but it 
does not appear that this has been the case. 

Carbon monoxide emissions 


Although the model does not currently predict stove emissions, 
some initial data have been sampled. Modified combustion 
efficiency is a commonly used metric for biomass cooking stove 
CO emissions. The calculation is given below in Eq. (19), where CO 
and C0 2 values are in mole/volume fraction. For the ideal goal of 
complete CO oxidation, combustion efficiency would be 100%. 
Modified combustion efficiency (MCE) averages for the 15 minute 
samples are shown in Fig. 11 with model predicted excess air ratio 
behavior. Measured excess air ratio behavior is shown with 
comparison to model predicted values in Fig. 10. 


MCE 


C0 2 

CO + C0 2 


(19) 


The initial results suggest that the excess air ratio is a promising 
tool for reducing carbon monoxide emissions. Higher combustion 
efficiency is observed at higher firepowers with correspondingly 
lower excess air ratios. Improved CO oxidation is likely explained by 
the reduced quenching and increased temperature at these decreased 
excess air ratios. The modified combustion efficiency trend is 
consistent between the two geometries, but offset to slightly higher 
firepowers for the larger stove. Excess air ratio is correspondingly 
offset, in agreement with the previous argument. 



Fig. 12. PM10 emissions factor for validation samples. 


Particulate matter emissions 

Particulate matter below approximately 10/ig for the 15-minute 
validation samples are shown in Fig. 12 in the form of an emissions 
factor ( EF ). EF is a way of normalizing PM emissions to firepower and 
is calculated using Eq. (20) shown below. 

_ g of Particulate Matter 
kg of Fuel Burned 

The initial data shown in Fig. 12 are interesting but more data is 
needed to make conclusions. All test runs exhibit lowered emissions 
factors in some medium firepower range; however, this range is not 
consistent between the four test runs, or within the two runs for each 
stove size. It is unclear whether the parameters of the model 
presented herein will be useful in understanding/improving partic¬ 
ulate matter emissions. 

An attempt to observe higher firepower/low excess air ratio behavior 

An upper limit to stove firepower is presumed in Section Stove 
firepower range to exist approaching low excess air ratios when, after 
accounting for poor mixing, as combustion becomes oxygen starved. 
Inefficient combustion is also expected in this oxygen starved region, 
providing an emissions upper limit in complement to the inefficient 
combustion region observed at low firepower and high excess air 
ratios in Section Carbon monoxide emissions. 

An attempt has been made to gather data in this region in order to 
observe these upper limits. As with previous validation up to this 
point, a cook piece was not in place. However, without a cook piece in 
place, it was observed that as the chimney region becomes oxidizer 
starved, combustion is able to proceed with entrainment air after 
exiting the chimney—the combustion chamber essentially becomes a 
torch. Operating in this mode, no detriment to CO modified 
combustion efficiency was observed. No attempt was made to sample 
PM in this case. Operation in this mode is not realistic for typical stove 
use and is not applicable as combustion proceeding after the chimney 
exit is not useful for cooking. 

However, this will not be the case in Part 2, Agenbroad et al. 
(2011), when sampling with the cook piece in place—practically 
applicable limits and interesting behavior will be observed operating 
in this region. 

Conclusions and future work 

The simple model developed herein has shown itself to be both 
accurate and practically useful in predicting bulk flow rate, temper¬ 
ature, and excess air ratio for a natural convection driven biomass 
cooking stove from geometry and firepower inputs—the simplified 
model predicts a sort of steady state behavior, and these steady state 
predictions have been shown to be in agreement with temporally 
averaged validation data sampled from typical stove operation. These 
stove flow prediction are intended as inputs for future model 
extensions predicting and improving cook stove heat transfer and 
emissions. Initial results suggest excess air ratio to be a useful tool for 
predicting and understanding carbon monoxide emissions as well as a 
stove’s firepower range. 

Part 2, Agenbroad et al. (2011 ), will expand on these results in the 
following ways: First, the two equation system summarized in 
Section The resulting two equation system will be further developed 
into a dimensionless form. Dimensionless form advantages include 
scale similarity and a reduction of experimentally independent 
parameters. Second, experimental validation will be extended to the 
more practical case of stove operation with a cook piece in place. An 
increased loss coefficient is experimentally determined, and the 
simplified stove flow model is shown to remain both applicable and 
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accurate. With the cook piece in place, high firepower/low excess air 
ratio behavior is both applicable and interesting. With the dimen¬ 
sionless form, Parts 1 and 2 emissions results combined, resulting in 
both expanded and increasingly convincing trends for both carbon 
monoxide and particulate matter. 

A discussion of future work in model validation as well as potential 
model applications to heat transfer and emissions will be deferred to 
Part 2, Agenbroad et al. (2011), conclusions. 

Acknowledgments 

This work has been supported by the Colorado State University 
Engines and Energy Conversion Laboratory with funding from 
Envirofit International and the Shell Foundation's Breathing Space 
Program. Special thanks to Melanie Sloan, Christian L'Orange, Dan 
Lionberg, and Dr. Anthony Marchese—for your help, expertise, and 
advice, but also for a great experience in general. 

References 

Agenbroad J, DeFoort M, Kirkpatrick A, Kreutzer C. A simplified model for 
understanding natural convection driven biomass cooking stoves—Part 2: With 
cook piece operation and the dimensionless form. Energy for Sustainable 
Developement 2011 ;15(2):169—75. 

Baldwin, S., 1987. Biomass stoves: Engineering design, development and dissemination. 
Ph.D. thesis, Center For Energy and Environmental Studies Princeton University, 
Princeton, New Jersey 08544 USA. 

Baldwin S, Geller H, Dutt G, Ravindranath NH. Improved woodbunning cookstoves: 

signs of success. Ambio 1985; 14(4/5 ):280-7 http://www.jstor.org/stable/4313163. 
Barnes DF, Openshaw K, Smith KR, Plas Rvd. The design and diffusion of improved 
cooking stoves. World Bank Res Obs July 1993;8(2): 119—41 http://www.jstor.org/ 
stable/3986529. 

Bruce N, Perez-Padilla R, Albalak R. Indoor air pollution in developing countries: 
a major environmental and public health challenge. Bull World Health Organ, 78(9). 
Tech. Rep.; 2000 http://www.sdelosp.org/sdelo.php?script=sci_arttext&pid= 
S004296862000000900004&lng=en. 

Bryden M, Still D, Scott P, Hoffa G, Ogle D, Bailis R, et al. Design principles for wood 
burning cook stoves. Tech. rep. 80574 Hazelton Road, Cottage Grove, Oregon 
97424: Aprovecho Research Center; 2005. 


de Koning HW, Smith KR, Last JM. Biomass fuel combustion and health. Bull World 
Health Organ 1985;63(1):11—26 http://www.ncbi.nlm.nih.gov/pmc/articles/ 
PMC2536350/pdf/bullwho00084-0024.pdf. 

Haines A, McMichael AJ, Smith KR, Roberts I, Woodcock J, Markandya A, et al. Public health 
benefits of strategies to reduce greenhouse-gas emissions: overview and implications 
for policy makers. Lancet 2010;374(9707):2104-14 http://www.sciencedirect.com/ 
science/article/B6TlB-4XST PDR-B/2/a228dee2ef027d9ad940cb5925b55el2. 

Heitor M, Moreira A. Thermocouples and sample probes for combustion studies. Prog 
Energy Combust Sci 1993;19(3):259-78 http://www.sciencedirect.com/science/ 
article/B6V3W-497B CKX-94/2/9358c4b23ald9a86f0d4677db577ce07. 

Hudelson, N. A., Bryden, K., Still, D., No Date. Global modeling and testing of rocket stove 
operating variations. Tech, rep., Aprovecho Research Center and Department of 
Mechanical Engineering, Iowa State University, 80574 Hazelton Road, Cottage 
Grove, Oregon 97424. 

Jetter JJ, Kariher P. Solid-fuel household cook stoves: characterization of performance 
and emissions. Biomass Bioenergy 2009;33(2):294-305 http://www.sciencedirect. 
com/science/article/B6V22-4T8H 38B-l/2/6e5ec727a3ce478944ab71c5010el409. 

Kee, R., Rupley, F., Miller, J., 1990. The chemkin thermodynamic data base. Technical 
Report SAND-87-8215B. Sandia National Labs., Livermore, CA(USA) (March). 

MacCarty N, Still D, Ogle D. Fuel use and emissions performance of fifty cooking stoves 
in the laboratory and related benchmarks of performance. Energy Sustain Dev 
2010;14(3):161—71 http://www.sciencedirect.com/science/article/B94T4- 
50RV8T5-l/2/4314261dcb920ad56be4bf0f65041f4b. 

Mwampamba TH. Has the woodfuel crisis returned? Urban charcoal consumption in 
Tanzania and its implications to present and future forest availability. Energy Policy 
2007;35(8):4221-34 http://www.sciencedirect.com/science/article/B6V2W-4NH6 
NFF-l/2/e3248896c6aal8a33cadd252d8b640dd. 

Naeher LP, Brauer M, Lipsett M, Zelikoff JT, Simpson CD, Smith JQKKR. Woodsmoke 
health effects: a review. Inhal Toxicol 2007;19(1):67-106 http://ehs.sph.berkeley. 
edu/krsmith/publications/2006%20pubs/JIT%20Woodsmoke2.pdf. 

Scott P. Updated basic stove design. Web; 2003 http://stoves.bioenergylists.org/ 
stovesdoc/Scott/rocket/updated%20basic%20stove%20design.pdf. 

Smith KR. Dialectics of improved stoves. Econ Polit Wkly 1989;24(10):517-22http:// 
www.jstor.org/stable/4394497. 

Smith KR. Health, energy, and greenhouse-gas impacts of biomass combustion in 
household stoves. Energy Sustain Dev 1994;l(4):23-9 http://www.sciencedirect. 
com/science/article/B94T4-4VC2 G5V-9/2/630493564774067b5b408ad2705b437a. 

Smith KR, Mehta S, Maeusezahl-Feuz M. The burden of disease from indoor air pollution 
in developing countries: comparison of estimates. Int J Hyg Environ Health 
2003;206(4-5):279-89 http://www.sciencedirect.com/science/article/B7GVY- 
4DS7 5H1 -7C/2/0523b72507d04000c3ba0b4ae602df3e. 

Tillman D. Wood as an energy resource. New York, NY: Academic Press Inc. January. 

Venkataraman C, Habib G, Eiguren-Fernandez A, Miguel AH, Friedlander SK. Residential 
biofuels in South Asia: carbonaceous aerosol emissions and climate impacts. 
Science March 2005;307(5714):1454-6 http://www.sciencemag.org/cgi/content/ 
abstract/307/5714 /1454. 


